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Abstract 


^ The  transient  pressure  on  a spherical  array  arising  from 
the  excitation  of  a simple  circular  element  has  been  obtained. 

Hie  cases  of  a Heaviside  step  and  a continuous  wave  (CW)  pulse 
were  considered.  The  low  frequency,  i.e.,  long  time,  response 
was  obtained  via  a classical  wave  harmonic  analysis.  High  fre- 
quency asymptotic  expansions  and  the  Watson  transformation  were 
used  to  isolate  the  short  time  transient  effects  into  efficient 
"creeping  wave"  series.  Dispersion,  caused  by  the  spherical 
curvature,  is  shown  to  produce  attenuation  of  the  traveling  wave 
and  theoretically  an  exponentially  decaying  pulse  of  infinite 
duration.  These  creeping  wave  series  are  not  valid  in  the  imme- 
diate vicinity  of,  and  on,  the  active  element.  However,  in  this 
region  the  Kirchhoff  approximation  has  been  shown  to  be  applicable 
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List  of  Symbols 


radius  of  sphere 


radius  of  piston  in  plane  baffle 

U/3 


I /2 


sound  velocity  in  acoustic  medium 

Green's  function 

Unit  or  Heaviside  step  function 

spherical  hankel  function  of  the  first  kind  of  order  m 


derivative  of  spherical  hankel  function  of  the  first  kind  of 
order  m 


u>/c 


u /c 
o 


Legendre  polynomial 


acoustic  pressure 
spherical  coordinates 

radius  vector  measured  from  center  of  piston  in  plane  baffle 

Laplace  transform  variable 

time 

time  of  maximum  pressure 


velocity 

density  of  acoustic  medium 

one-half  angle  subtended  by  piston  in  spherical  baffle 


root  of  the  derivative  of  the  Airy  function 


circular  frequency  of  CW  pulse 


indicates  Fourier  or  Laplace  transform 


1.  Introduction 


The  spatial  fluctuations  of  surface  pressure  on  a large  array  give  rise 
to  a surface  pressure  distribution  whose  average  over  individual  transducer 
elements  varies  widely  from  one  transducer  element  to  smother,  if  these  ele- 
ments are  small  in  terms  of  wavelengths.  Consequently,  unless  the  differences 
in  radiation  loading  are  electrically  or  mechanically  compensated,  the  dynamic 
configuration  of  the  array  and  hence  the  sound  field  will  not  be  as  specified 
by  the  beam-forming  calculations.  This  difficulty  does  not  arise  if  elements 
have  characteristic  dimensions  of  one-half  wavelength  or  more  because  the  fluc- 
tuations in  surface  pressure  amplitude  tend  to  average  out  when  integrated  over 
the  surface  of  an  element.  The  need  of  compensating  for  surface  pressure  fluc- 
tuations is  therefore  primarily  associated  with  the  long  wavelengths  resulting 


from  the  recent  trend  toward  long-range  and  hence  low-frequency  echo-ranging. 
The  need  and  implementation  of  compensating  for  the  variation  in  the  radi- 


ation loading  has  been  thoroughly  analyzed  . For  this  purpose  the  radiation 
impedance  matrix,  i.e.  the  surface  pressure  averaged  over  individual  transducer 


elements  when  a single  transducer  element  is  active,  must  be  computed. 

This  study  provides  an  analytical  solution  suitable  for  computing  the 
transient  radiation  loading  on  a spherical  BQS-6  type  array  at  the  onset  of 
a step  or  a continuous-wave  (CW)  pulse,  with  no  restriction  as  to  beam  steering 
and  formation.  Preliminary  numerical  results  were  presented  at  the  Fall  19<>9 
meeting  of  the  Acoustical  Society  of  America  in  San  Diego  (paper  UBl). 


The  required  surface  pressure  can  be  analyzed  by  means  of  the  classical 

2 

wave  harmonic  formulation  . However  this  series  exhibits  poor  convergence 
for  large  arrays  in  terms  of  wave  length,  particularly  in  the  vicinity  of  the 
transient  vave  fronts. 


An  alternative  approach,  which  is  undertaken  in  this  work,  had  origi- 

3 

nally  been  developed  by  Watson  to  compute  the  diffraction,  around  the  earth 
of  electromagnetic  waves  from  a point  source.  That  formulation  has  re- 
cently been  applied  to  the  evaluation  of  the  surface  pressure  on  large 
4 5 

sonar  arrays  * . This  "creeping-wave"  formulation  is  tantamount  to  an  ex- 
pansion of  the  pressure  in  negative  powers  of  ka  and  is  therefore  partic- 
ularly well  suited  to  large  arrays.  Unfortunately,  it  does  not  yield  the 
self- impedance , i.e.,  the  pressure  on  the  active  element.  However,  in  the 
large-ka  range  where  Watson's  technique  is  useful,  the  self- impedance  can 
be  approximated  by  using  the  self-impedance  of  a plane  piston  in  a plane 
baffle,  with  a minor  correction  applied  for  baffle  curvature  to  the  resis- 
tive component.  For  the  ka  and  transducer  element  dimension  of  the  BQS-6, 
the  resulting  correction  to  the  self- impedance  is  shown  to  be  negligible. 

2.  Mathematical  Model 

Consider  the  transient  radiation  loading  associated  with  a single, 
active  element  located  on  a spherical  array  when: 

(A)  the  active  element  undergoes  a velocity  step 

w(t)  = WH(t)  (1] 

(B)  the  active  element  initiates  a "CW-pulse" 

w(t)  = WH(t)  sin  (i)Qt  (2] 

The  transient  surface  pressures  thus  constitute  the  transient  equivalent  of 

the  familiar  radiation  impedance  matrix.  Item  (A)  constitutes  an  extension 

6 

to  a spherical  baffle,  of  Miles  analysis  of  a piston  in  a plane  baffle  . 
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Item  (B)  constitutes  a similar  extension  of  Mangulis'  analysis  . 

For  the  BQS-6  the  array  is  large  in  terms  of  the  CW-wavelengths  (e.g. 
ka  ■ 32),  and  the  active  element  is  small  (e.g.  a linear  dimension  of  lU\/l*5, 
corresponding  to  a subtended  angle  of  3.5°). 

Under  these  conditions,  the  surface  pressures  can  be  computed  over  mof.t 
of  the  sphere  as  though  the  active  element  were  a point  source . This  assump- 
tion is  not  valid  in  the  immediate  vicinity  of,  and  on  the  active  element. 

The  spherical  polar  axis  (0  = 0)  is  made  to  coincide  with  the  center  of  the 
active  element  (Figure  l).  The  active  element  subtends  the  angle 


0 < 0 < a 


The  mathematical  problem  may  therefore  be  stated  as  the  solution  to  the 
wave  equation  in  spherical  coordinates 


2 2 
VsPh  ? + k P * 0 


(3) 


subject  to  the  boundary  condition 


w (a,0;t)  = f (t)  H (a  - 0) 


(10 


where  f (t)  is  defined  by  equation  (l)  or  (2)  and  k = w/c 
3.  Wave- Harmonic  Analysis 

The  transient  radiation  loading  can  be  formulated  in  terms  of  the  inverse 
Fourier  transform  of  the  transient  velocity  distribution 

h (ka) 

p (a,6;t)  = l An  Pn  (cos  0)  / ir  («>)  ^TkaT  e 
n*0  — 00  n 


-iwt 


du 


(5) 


v.a! 


________ 
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Expanding  this  velocity  distribution  in  Legendre  functions0  the  prescribed 
boundary  condition  becomes 

00 

w (a,0;t)  * f (t)  y [P  . (cos  a)  - P ..  (cos  a)]  P (cos  0)  (6) 

n-1  n+l  n 

n=0 

where  f (t)  is  defined  by  equation  (l)  or  (2).  If  the  piston  radius  is • 


allowed  to  tend  to  zero. 


w (0  ;t  ) 3 \ a2  f (t)  l (n  + h P (cos  0 ),  a2  « n~2 


o*  o'  2 


2'  n 


For  any  finite  piston,  however  small,  the  higher  order  terms  must  eventually 
be  written  in  the  form  of  Equation  6,  when  n becomes  large  enough.  This  un- 


fortunately, will  be  shown  to  be  insufficient  to  insure  convergence  of  the 
surface  pressure  at  the  beginning  of  the  pulse. 


The  Fourier  transform  of  the  pressure  therefore,  is  a series  in  the  form 


$(a,0;u>)  = - f («)  l [Pn_1(cos  a)-Pn+1(cos  a)]  ^'(k'a)  pn^cos  ^ 

n=0  n 


h (ka) 


The  corresponding  inverse  transform  can  be  evaluated  using  residue  theory 
which  expresses  the  pressure  as  series  of  wave-harmonics.  The  transient  com- 
ponent of  the  solution  is  associated  with  the  residues  arising  from  the  roots 

of  h'(ka)  = 0.  The  number  of  roots^  and  hence  of  residues  associated  with 
n 

the  n^*1  wave-harmonic  is  (n+l).  Any  steady-state  contribution  is  produced 
by  the  singularities  of  f (to).  The  transient  pressure  field  thus  becomes: 

For  case  (A) 

n h (z  )exp(-iz  ct/a) 

p(a,0;t)  = ipctfH(t)  J £ ' z ') X 

n=0  J=1  nj  n nj ’ 


IP  , ( cos  a)  - P .. (cos  o)]P  (cos  0) 
n— 1 n+i  n 


I 


♦ * -V 


— 
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For  case  (B) 


ln+1  k a h (z  )e_iznjct^a  ..  h (k  a)  » 

p(a,0;t)=pc*H(t)  l < [ — -n  ^ ^ + 2Tk~aT  h^(k°a)  e °*  a 

n=0  J=1  z2  -(k  a)2  hn  Unj'  2Uoa;  Vo' 


,1 


nj  ' o 


[P  , (cos  a)  - P ..(cos  a)]  P (cos  8) 
n-1  n+1  n 


(10) 


where  zqj  are  the  complex  roots  of  the  equation  h/(ka)  = 0 and  kQ  = u)q/c. 

The  first  term  of  equation  (10)  represents  the  transient  effect.  The  second 

term  arises  from  a pole  at  to  = u and  is  the  steady-state  contribution. 

o 

Consider  the  convergence  properties  of  the  residue  series,  Equation  9 and 
10.  The  surface  pressure  is,  by  virtue  of  the  Initial  Value  Theorem 


p(a,0;t)  = 
lim  t-*-  0 


- itop(a,0 ;<o) 
lim  to  +■  00 


(11) 


Letting  k -*■  00  in  Equation  8,  the  ratios  h^/h^  tend  to  i.  Substituting  in 
the  above  asymptotic  relation,  the  early-time  surface  pressure  becomes  pro- 
portional to 


p(a,6;t)  £ [P  (cos  a)  - P (cos  a)]  P (cos  0) 

. . ...  __  n-1  n+1  n 

lim  t ->-0  n=0 


(12) 


This  series  converges  no  faster  than  the  series  representing  the  velocity 
distribution  Equations  6 and  7.  The  wave-harmonic  series  is  therefore  not 
a practical  formulation  for  evaluating  the  pulse  at  early  times  particularly 
for  small  pistons,  for  which  Equation  7 applies  to  a large  number  of  terms. 
It  thus  becomes  apparent  that  the  early-time  transient  radiation  loading  on 
a piston  set  in  a spherical  baffle  is  not  tractable  by  means  of  the  rigorous 
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wave-harmonic  formulation  described  above.  This  is  in  fact  true  for  any  values 


of  the  (0,t)  domain  corresponding  to  the  vicinity  of  the  propagating  transient 
wave  front.  This  is  in  contrast  to  the  circular  piston,  whose  double-transform 
• solution  yields  the  radiation  loading  at  early  times  without  difficulty. 

Asymptotic  techniques  are  therefore  required.  In  the  next  sections 'an 
approximate  solution  is  developed  for  obtaining  the  transient  pressure  on  the 
array  specialized  to  the  situation  of  slowest  convergence,  i.e.,  the  limiting 
case  of  the  small  piston.  Furthermore,  the  solution  is  deliberately  special- 
ized  to  early  times,  i.e.,  to  the  high-frequency  end  of  the  spectrum.  This 
is  achieved  by  applying  the  Watson  transformation  to  the  pressure  transform, 
in  Equation  8,  specialized  to  small  a. 

U.  Creeping-Wave  Formulation  - Velocity  Step 

The  following  analysis  uses  Watson's  transformation  to  replace  the  steady- 
state  wave-harmonic  series  with  a residue  or  "creeping-wave"  series  which 
captures  the  high-frequency  components  of  the  solution  in  its  leading  terms. 


The  resulting  transform  takes  a form  suitable  for  integration  by  steepest- 
descent  approximation.  One  thus  achieves  not  only  rapid  convergence,  but  one 
has  cast  the  integrand  in  a form  which  permits  analytical  evaluation  of  the 
inverse  Fourier  transform. 

This  technique  has  been  validated  for  the  CW  surface  pressure  generated 
on  a large  cylinder  by  a point  source"^.  In  this  case,  the  pressure  is  a 
Fourier  transform  over  all  axial  wavenumbers,  rather  than  over  frequency, 
but  the  same  asymptotic  integration  applies. 


Consider  case  (A),  i.e.,  the  velocity  step.  The  formal  expression  for 
the  transient  surface  pressure  specialized  to  infinitesimal  pistons  is  ob- 
tained by  using  Equation  7,  in  lieu  of  Equation  6,  for  the  modal  velocity 


— 


coefficients.  Thus  using  the  small  piston  approximation  to  expand  the  spatial 
dependence  of  Equation  4,  i.e.  Equation  7,  and  the  Fourier  transform  of  Equation 


1,  Equation  5 becomes: 

p(a,0 ;t ) = - pC  W / 


2 /.  00  -iwt 
e 


, h„(ka) 

l <n  * 5>  i%?  pn<‘°8  6)  4“ 

n=0  n 


(13) 


Note  that  the  integrand  of  Equation  13  apparently  contains  a singularity  at  w = <>• 
However,  since 


h (x) 


Limit  , ,r 

*K>  V*1 


= 0(x) 


the  integrand  is  finite  at  w = o and  in  fact  no  singularity  exists.  The  wave- 
harmonic  series  in  the  above  expression  can  be  written  as  Watson's  residue  series: 


n=0 


h (x) 

: V7 — r P (cos  0)  = 
h (x)  n 
n 

Y 1 

iv  0 
e m 

Jo  K71 

(sin  0)' 

: - 21/3, 

* 

x5'3  exp 

1/2 


7/6 


iv  (2tt-6) 
e m 


(i4) 


(sin  0)^^  [sin(2ir  0)]"^2 


m=0 


iv  IT 

e m 


m' 


(15) 


where 


x = ka 

v = x + |a'|  (-x/2)1/3 
m ' m* 

- |a^|  = roots  of  the  derivative  of  the  Airy  function^ 

= 1.02,  3.25,  4.82,  ....  [3ir(4m  + l)/8]2/3 
The  rapid  convergence  of  the  creeping-wave  series  for  large  x,  i.e.,  at  high 
dimensionless  frequencies  ka,  is  apparent  if  the  exponent  is  split  into  its 
real  and  imaginary  components,  thus  displaying  an  exponentially  decaying 


(16) 


factor : 


iv  0 

e n ' exp 


1/3 


I 


f 


1 

I 


It  is  precisely  this  good  convergence  of  the  series  at  higher  frequencies 
vhich  makes  it  suitable  for  exploring  the  early  times  of  a transient  event. 

The  creeping-wave  series  does  not  converge  on,  or  in  the  immediate  Vicin- 
ity of,  the  active  element.  Fortunately,  in  the  high-frequency  limit  where 
this  formulation  is  useful,  the  surface  pressure  in  this  region  can  be  approx- 
imated by  ignoring  baffle  curvature  as  0 a.  The  steady-state  surface  pres- 
sure field  generated  by  a piston  on  a plane  baffle  blends  smoothly  with  the 
field  computed  from  the  creeping-wave  series^.  Similarly,  under  transient 
conditions , the  surface  pressure  on  the  active  element  will  be  shown  to  be 
approximated  at  early  times  by  the  solution  derived  for  the  plane  piston 
(section  7). 

In  the  early-time,  high-frequency  limit,  the  (2ir  - 6)  term  in  Equation 
lU  can  be  dropped,  except  in  the  immediate  velocity  of  the  antipodes  0 = ir. 
Thus , retaining  only  the  0 - term  and  substituting  the  resulting  creeping-wave 
series  in  Equation  13,  the  transient  surface  pressure  becomes 


„ ...  2 7/6  -i*/  12 

, n \ p Wa  a e 

p(a,0;t)  = i T—rrzr 

8(2*  sin  0p  c1/b 


\ exp  (-iwt  + iv  0)  die,  0 ^ 0,  * 

m=0  Cm  -°°  171 


(17) 


Since  which  is  of  order  wa/c,  is  large  in  the  high-frequency  limit 
the  phase  angle  is  caused  to  vary  rapidly  with  u.  The  inverse  transform  in 
Equation  17  thus  lends  itself  to  integration  by  the  method  of  steepest  decent 
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whereby  an  integral  of  the  form 

Is/"*  (»)  ei’J'(u)  a,,  (18) 

—CO 

can  be  approximated  by  the  contribution  of  the  saddle-point  to  defined  by  the 
condition 

= 0,  to  = to  (19) 


This  yields 

I - (2*0**  »(5)  ei^,(u>) 


(20) 


Using  Equation  IT,  Equation  19  yields 


d\|t 


3v 


— = - t + 6 t*—  = 0,  to  = to 

a to  a to  ’ m 


(21) 


From  Equation  16 


8v  m a |fm|  (- a )^3 

3“  =c  3 (2c«o2)"1/3 


(22) 


When  this  is  substituted  in  Equation  21,  one  can  solve  for  the  saddle-point 


value , to  = to 


m 


(o'  6 /3)3/2  v 

to  = o/p  (-a/2cr 

lt-(a0/c)]3/2 


(23) 


As  [t-(a  0/c) }*0,  to  ■*  ».  The  condition  t •+  (a  6/c)  therefore  represents  a 


m 


wave  front  spreading  over  the  sphere  from  0 = 0 , by  virtue  of  the  Initial 


Value  Theorem.  The  stationary -phase  approximation  enveloping  the  spherical 


'baffle  Integral  oust  therefore  be  multiplied  by  the  coefficient  H [t-(a  6/c)] 


The  modulus  of  the  Integrand  in  Equation  17  becomes 


( |<*lI  6/3) 


l/u 


♦(«  ) 

" [t-(a  6/c)] 


, >1/12  is  /12 

37^  ( a/2c ) e 


(24) 


The  second  derivative  of  the  phase  ang  is  also  required  in  constructing 
Equation  20 


a2  * 


m 


2 

9 v 


6 


3u> 


3w 

,2/3 


a'  0 
m _ 


L 


(25) 


- .»  - » 


a17^  ( 


m 


e) 


m 


Finally,  the  phase  angle  evaluated  at  to  = is  required: 


ij;  (to)  = 
m 


- <u(t-  — ) + 
c 


m’ 


. ua>l/3 
2c' 


i(2a)1/2(|aj;|0)3/2 
33/2  cl/2[t_(ae/c)]l/2 


, u>  = 


0) 

m 


(26) 


The  phase  angle  is  purely  imaginary,  thus  giving  rise  to  an  exponentially 
decaying  term.  When  the  results  in  Equation  24,  25,  and  2 6 are  substituted 
in  Equation  20,  the  stationary-phase  approximation  to  the  integrals  in 
Equation  17  is  obtained: 


/ 


dui  = 


_1/2  . 

IT  O 


eal/3  ,-i.  /12 


(2«)1/2( 


1 /o  0/9  < 

1/3[t-(a0/c)]3/2 


m 


e ) 


3/2 


33/2c1/2[t-(.e/c>]1/2 


I 


H(t-  Si) 

C 


(27) 


When  this  result  is  substituted  in  Equation  17,  and  expressing  |a^|  in  terms 

J.  /o 

of  the  coefficients  c^  = |a^|/2  ' , the  expression  for  the  transient  pressure 


finally  becomes: 


p(a,0;t) 


pW  q2a3/^2e  H[t-(a9/c)1  \ 

4 (6c  sin  e)1^2  [t-(a0/c ) ]3/2  m=0 


U(2a)1/2  (cm0)3/2 
33/2c1/2  [t_(a0/c)]l/2 


This  solution  is  not  applicable  at  or  near  0=0.  Here  it  may  supplement 
the  solution  for  the  impulsively  accelerated  plane  piston  derived  by 
Miles**  and  reproduced  in  section  6. 

It  was  mentioned  earlier  that  the  creeping-wave  solution  does  not  con- 
verge in  the  illuminated  region  a > 0.  Actually,  the  validity  of  the  solu- 
tion is  more  generally  limited  in  the  region  near  0=0  and  0 = it  because 
it  embodies  a large  |v  sin  0|  approximation  to  the  hypergeometric  function 

P , /o,  a function  which  arises  from  P as  a result  of  the  Watson  trans- 
v-  1/2’  n 

formation.  Consequently,  the  creeping-wave  solution  predicts,  incorrectly, 

a surface  pressure  larger  than  the  plane  baffle  solution  in  a region  bounded 

2/5 

approximately  by  0 < 50  o , with  both©  and  a in  degrees.  For  a comparable 
active  element,  this  region  is  far  wider  than  for  the  steady-state  situation. 
In  the  transient  situation,  the  a-dependence  of  the  extent  of  this  region 
arises  because  the  creeping-wave  solution  which  predicts  a pressure  varying 

• . * 2i 

linearly  with  the  volume  velocity  Q (=  Wo  } , is  smaller  than  the  small-piston 
(or  large-range)  limit  for  the  plane  baffle,  which  varies  as  Wb,  i.e.  as  Q/b, 


and  therefore  becomes  singular  as  h •*  0.  This  difference  arises  because  for 
field  points  on  a curved  baffle,  the  acoustic  shadow  cast  by  the  baffle  com- 
pounds the  spherical  spreading  loss  thus  eliminating  the  singularity  associ- 
ated with  a point  source.  For  a plane  baffle,  where  the  shadow  effect  is 

absent  spherical  spreading  alone  mitigates  but  does  not  eliminate  this  singu- 

• * 2 

larity.  On  the  active  element,  a singularity  of  the  form  Q/h  arises. 

Another  significant  difference  between  the  two  baffle  configurations  arises 
from  the  dispersive  and  attenuated  character  of  the  creeping-waves:  For  plane 
baffles  the  pressure  pulse  displays,  like  the  piston  velocity,  a time-deriva- 
tive which  is  infinite  at  the  leading  edge.  For  the  curved  baffle,  the  onset 
of  the  pulse  represented  by  Equation  28  is  gradual,  all  time  derivatives  being 
zero  near  the  leading  edge,  as  expected  from  the  fact  that  exponential  attenua- 
tion increases  monotonically  with  frequency  thus  extinguishing  the  signal 
associated  with  the  leading  edge  of  the  pulse.  Furthermore , while  the  pulse 
duration  T = 2b /c  tends  to  zero  for  the  plane  baffle  as  the  source  dimension 
2b  tends  to  zero,  this  is  not  the  case  for  the  curved  baffle.  The  reason  is 
the  dispersive  character  of  the  creeping-wave  which  lengthens  the  pulse  dura- 
tion as  the  pulse  envelops  the  sphere.  The  time  at  which  the  maximum  of  a 
given  creeping-wave  mode  occurs,  obtained  by  setting  the  time  derivative  of 
Equation  28  equal  to  zero,  is 


Since,  for  all  cases  of  interest,  only  the  lowest  creeping-wave  mode  is 
required  at  the  pressure  peak,  one  obtains  an  analytical  expression  for 
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- I ' 


the  pressure  maximum: 


p(a,0;tQ) 

• 2 
p cW  a 


37  e"3 


Fla' 


9/2(6T  sin  9)1/2 


6.20 

( 0T  sin  e)1^2 


(30) 


Thus,  the  further  the  angular  distance  traveled  by  the  pulse , the  more  the 
peak  pressure  falls  behind  the  wavefront.  Furthermore  until  the  caustic 
at  0 — it  is  approached  the  spreading  loss  far  exceeds  the  spherical  spreading 
loss  observed  on  the  plane  baffle.  A final  comment  on  the  solution  is  that, 
even  though  a high-frequency  asymptotic  technique  has  been  used,  the  solution 
tends  to  the  correct  late-time  limit,  viz.  zero. 

The  result  in  Equation  28  coincides  with  the  solution  derived  by 
Friedlander11  using  a different  formulation.  His  Green's  Function  approach 
is  described  in  the  next  section  and  used  to  obtain  the  transient  pressures 
at  the  beginning  of  a CW-pulse.  Even  though  Friedlander' s asymptotic  approx- 
imations and  method  of  integration  are  identical  to  those  above,  his  procedure 
differs  in  that  he  does  not  apply  Watson's  transformation  to  the  wave-harmonic 
series , but  finds  a solution  in  the  form  of  a creeping-wave  series  from  the 
wave  equation. 

5.  Green's  Function  Formulation  - CW-Pulse 

The  transient  pressure  expressed  as  an  inverse  Fourier  transform  (Equation 
5)  can  be  written  in  terms  of  the  steady-state  Green's  function,  specialized 
to  surface  pressures  (R  = a),  and  axisymmetric  velocity  distributions.  This 
Green's  function  will  be  designated  as  G,  to  distinguish  it  from  the  transient 
or  impulse  Green's  function  G used  by  Friedlander.  The  steady-state  pressure 
will  be  similarly  expressed  as  f,  as  will  the  displacement  distribution  $(0;u) 
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of  the  boundary.  The  CW  surface  pressure  can  now  be  written  as 


2 2 t 

$(a,8;u)  ■ 2ir  pu  a / G(a,0  ;w)  tf(0  ;u)  sin  6 d6 

q o o o o 

(31) 

2 2 

s - ir  pa  a &(a,0|a,O;o>)  V(io),  a < < 1 

The  latter  approximation  is  not  necessary,  but  simplifies  the  analysis 
without  impairing  the  results  in  the  region  where  the  creeping-wave  solu- 
tion is  efficient. 

Noting  that  Equation  5 is  of  the  form 

p(a,0;t)  » ff  (u)  e“  du>  (32) 

it  (ai) 


and  expressing  the  ratio  f/it  in  terms  of  Equation  31,  the  transient  pressure 
becomes 


p(a,0;t) 


2 2 
pa  a 


/ G(a,0|a,O;w)  t)(u))  e“iaVt  da) 


(33) 


For  the  velocity  history  in  Equation  2;  i.e.,  the  CW-pulse, 


V(t)  = t u H (t)  cos  to  t 
o o 


(3U) 


whose  Fourier  transform  is 


V(w) 


i W b>  u 
0_ 

2 2 
u - 0) 

o 


(35) 


When  this  is  substituted  in  Equation  33,  the  transient  pressure  becomes 


p(a,6;t) 


iW  to  / 


u 


G(a,0[a,O;u)  e”  u du 


(3< 


u - to 


The  resulting  pressure  in  this  situation  is  made  up  of  two  contributions  to 
the  inverse  transform.  A saddle-point  contribution  which  represents  the 
transient  effect  will  be  shown  to  tend  to  zero  with  increasing  time.  The 
residue  at  u>  ■ <*>o  is  not  attenuated  with  time  and  in  fact  is  the  steady- 

state  creeping-wave  solution. 

A. ) Steady-State  Response 

The  steady-state  problem  of  a point  source  on  a spherical  radiator  has 
previously  been  solved"^.  Fbr  the  sake  of  continuity,  the  results  are  pre- 
sented below 


p(a,e)  = P-  &(ka) 


13/6 


(32ir  sin  e)l/2a3 


I r 

m=0  m 


[-a  (ka)1/30+i4>  (0)] 

m m 


0 / 0,  ir 


(3' 


pc2Q  Py  . (-cos  e)(ka)®^3 

p(a,e)  « 

Ua3c 

o 

{-a  (ka)1/3*  -iut  +i*  [ka  + c (ka)1/2  + 2/3]} 

„ ° o . 

e 0 * IT 


(3( 


where 
Q * Vita2 

^ " 12^  ♦ ka  [l  ♦ cffl(ka)”2'3]  0 - ut 

B. } Transient  Response 

Consider  the  steepeat-decent  evaluation  of  the  integral  arising  from 


MM 


the  high  frequency  end  of  the  spectrum.  R>r  finite  the  factor 


does  not  affect  the  saddle-point  value  arising  from  the  asymptotic  evalua- 
tion of  6 (a,e| a,0;ai) . The  inverse  may  therefore  he  obtained  by  inverting 
the  transform  of  the  Green’s  Function  and  modifying  the  result  by  the  above 
multiplicative  factor  evaluated  at  the  resulting  saddle-point  value. 

The  integral 

QO 

f G (a,0(a,O;oj)  e dm  (3! 

—00 

is  by  definition  of  the  inverse  R>urier  transform,  2ir  times  the  impulse 
Green's  function  G (a,0 ,t [ a,0 ;0) . 

Friedlander  obtained  an  asymptotic  expression  for  this  Green's  function 
for  0 / 0,  ir  and  for  0 = it  using  one  creeping-wave  mode.  Using  consistent 
notation  and  expressing  the  Green's  function  in  terms  of  the  roots  ( ~ 
the  derivative  of  the  Airy  integral  these  expressions  are 


G(a,ir;t|a,0;0) 


air 


9 A 

-J  H (ct  - air)  X 


As  explained  earlier  a different  formulation  must  be  used  for  0 s 0.  To 
obtain  the  pressure  close  to  the  active  element,  Friedlander's  result  must 
be  modified  by  retaining  more  than  one  "creeping-vave"  mode.  The  number  of 
modes  required  as  a function  of  location  6 and  time  will  be  discussed  in  a 
later  section. 

Vhen  the  transient  pressure  is  formulated  as  an  inverse  Laplace  transform, 
the  procedure  consists  of  multiplying  the  terms  in  Equations  Uo  and  4l  by 
the  Laplace  transform  of  Equation  3^,  evaluated  at  the  saddle-point.  Friedlander 
who  formulates  his  impulse  Green's  function  as  an  inverse  Laplace  rather  than 
Fourier  transform,  gives  the  saddle-point  value  of  the  Laplace  transform  var- 


iable "s"  as 


/ a \l/2  [ “m8  1: 

2c  ([3  t-(a0/c)]J 


To  use  this  result,  ve  require  the  Laplace  transform  of  Equation  3^: 


Substituting  the  saddle-point  value  s = s , this  becomes 

m 


V (s  ) 21/2[3(ct  - a6)  “ 0]3/2  (k  a)2 

m m o 


Multiplying  the  terms  in  Equation  1*0  by  this  quantity,  the  transient  pressure 
becomes 


i 


p(a,0;t) 


3pcW*2  (k  a)205/2 


4(sin  0) 


•H  (t  - Si)  X 


j °.3/g  ^-^>1/2<V>3/2i 

m«l  5U(u  t - k afl  )3  + k a(o  0)3 


The  pressure  is  finite,  even  at  the  wavefront.  However,  at  the  center  of 
the  shadow  zone,  the  focusing  action  is  strong  enough  to  generate  a pressure 
singularity 

„lA_15A/,  n2  *_2  a 3/U  aiT\ ' 

3 n (k  a)  pcWa  —p-  H(t  - — ) 

O / . \J/4  C 

_/-  - j-x  _ (ct  - air) 

pla.ifjtj  o x 


m*l  [51*(u  t - k a6)3  + k a(a  0)3] 
o o o m 


6.  Plane  Baffle  Analysis 

As  previously  mentioned,  the  creeping-wave  formulation  does  not  yield 
the  pressure  on  the  active  element  and  is  inefficient  in  its  immediate 
vicinity.  In  this  and  the  proceeding  sections  the  plane  baffle  analyses 
of  Miles  and  Mangulis  are  presented  for  completeness,  and  shown  to  be  appli- 
cable on  and  near  the  piston. 

Miles^  obtained  the  solution  for  the  transient  pressure  on  a circular 
piston  embedded  in  an  infinite  plane  baffle.  The  exciting  pulse  was  assumed 
to  be  a Heaviside  step  function.  Transform  techniques  were  used  and  in  this 
case  a closed  expression  was  available  for  the  inverse  transform. 


I 

! 


H 1 


ct  < (a  - r) 


p(r ,0;t) 


2 2 A 2 v2 
c t + r - b 


it  cos 


(r  - b)  < ct  s (r  + b) 


(r  + b)  < ct  < ( r - b) 


where  r is  the  radial  distance  measured  on  the  plane  of  the  baffle  and  b 
is  the  piston  radius. 

7 

Mangulis  studied  the  case  of  a pulsed  CW  excitation  and  was  able  to 
obtain  an  expression  for  the  resulting  force  on  the  piston 


p(r,0;t)  = ^{We^o^tl-  — j sin^  0 


-i2  k b cos  © 


dG]>  0 < t < Q- 
~ c 


-1  ct 
cos  2b 


, . J,(2kb)  iH.  (2k  b) 

p(r,0;t)  = Im{We  Wo  [l j— g2 — + — k b ° — ]}  t > 2b/c 

o o 


where  J,  is  a Bessel  function,  H,  is  a Struve  function,  and  k = to  /c  . 

1 1 o o 

7.  Correction  for  Baffle  Curvature  on  and  Near  Piston 

The  effect  of  the  baffle  curvature  on  the  maximum  pressure  sensed  by 

the  piston  and  by  the  immediately  adjacent  baffle  area  elements  can  be 

shown  to  be  non-existent  for  the  former  and  small  for  the  latter.  Immediately 

after  piston  acceleration,  when  the  entire  piston  senses  the  maximum  pressure, 

pcW,  the  pressure  is  confined  to  the  piston.  Therefore  the  piston  "does  not 
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know"  whether  it  is  in  a plane  or  curved  baffle,  since  the  latter  has  not 
interacted  with  the  shock  wave.  The  immediately  adjoining  region  of  this 

baffle,  senses  its  maximum  pressure,  pcW/2,  immediately  after  piston  accel- 

* 

eration  , before  more  remote  regions  of  the  baffle  have  been  reached  by 
the  shock  wave.  This  pressure  value  should  be  affected  to  a negligible 
extent  by  baffle  curvature,  as  long  as  the  angle  between  the  plane  of  the 
piston  and  the  plane  tangent  to  the  baffle  area  element  is  extremely  small. 

It  is  only  as  the  shock  wave  reaches  more  remote  regions,  where  the 
shock  front  has  been  rotated  through  increasingly  large  angles , that  baffle 
curvature  may  be  of  significance.  It  will  now  be  shown  by  semi-quantitative 
reasoning  that  this  effect  is  negligible  in  the  entire  range  of  0 where  the 
plane  baffle  solution  rather  than  the  creeping-wave  solution  must  be  used. 

If  the  empirical  steady-state  correction  for  baffle  curvature  is  intro- 
duced in  the  inverse  Laplace  transform,  analytical  integration  becomes 
impossible.  The  size  of  the  correction  which  would  be  obtained  if  the 


•More  precisely  the  maximum  pressure  is 


max 


pcW  -1  r,_  a2, 1/2, 

~ cos  I C1  " T)  J 


pcW 


as  r + a 


pcW 

= £ — as  r >>  a . 
nr 

This  maximum  pressures  occurs  at  time 

. 1/2  2,1/2 

m c 

The  sudden  drop  from  a maximum  pressure  pcW  everywhere  along  the  piston  to 
half  that  value  on  the  Immediately  adjoining  baffle  assumes  no  viscosity, 
i.e.  an  acoustic  boundary  layer  of  zero  thickness. 
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baffle  curvature  was  taken  into  account  by  evaluating  the  curvature- 
corrected  inverse  transform  numerically  will  be  seen  to  be  negligible  from 
the  following  consideration 


The  pressures  on  the  spherical  baffle  and  on  the  plane  are  approximately 


related  as 


p s p -1 

sphere  plane  [l  - (ka)  ] 


The  relevant  value  of  w,  or  ka,  is  the  saddle-point  value  Equation  23, 

corresponding  to  the  time  of  maximum  pressure  (t  ),  i.e. 

o 


j c 


When  this  value  of  t is  substituted  in  Equation  23,  one  finds 


for  t = t 


with  0 in  radians.  This  quantity  is  so  large  within  the  range  of  0 where 
the  curvature-corrected  plane-baffle  solution  is  applicable  as  to  make  the 
correction  term.  Equation  1+9,  negligible 
8.  Computation 

The  region  of  the  (0,t)  domain  in  which  the  wave  harmonic  series 
(Equations  9 and  10 ) are  applicable  is  limited  by  the  difficulty  in  obtain- 
ing the  values  of  h"  (z)  evaluated  at  the  (n  + l)  complex  roots  of  the 
m 

equation  h'  (z)  = 0.  To  extend  the  number  of  roots  available  beyond  exis- 
ting tabulated  values  (Table  l)  may  require  the  use  of  a computer  subroutine 


which  obtains  the  complex  zeros  of  an  arbitrary  function.  A qualitative 


measure  of  the  number  of  terms  required  as  a function  of  position  and  time 
may  be  obtained  from  consideration  of  the  stationary  phase  value  of  ka  given 
by  Equation  23.  In  the  CW  case,  the  extent  of  applicability  of  the  series 


i 


t 


is  also  proportional  to  u>o. 

The  creeping-wave  series  (Equation  28  and  1*5)  for  the  transient  pressure 
exhibit  excellent  convergence  in  its  domain  of  applicability,  i.e.,  in  the 
vicinity  of  the  wave  fronts  away  from  the  active  element.  In  most  cases,  a 
maximum  of  four  terms  has  been  found  to  satisfy  a five  precent  convergence 
criterion.  For  0 » a only  the  first  creeping-wave  mode  is  required. 

Equations  28  and  1*5  may  also  be  used  to  evaluate  the  transient  pressure 
when  the  active  area  can  not  be  considered  a point  source.  This  situation 
would  exist  if  a combination  of  elements  were  excited  or  if  the  creeping-wave 
formulation  was  desired  in  the  immediate  vicinity  of  the  active  area.  In 
this  case.  Equations  28  and  1*5  may  be  considered  as  Green’s  functions  and 
the  resulting  pressure  obtained  via  a numerical  integration  over  the  active 
area.  This  is  described  in  Figure  2. 

9.  Summary 

The  transient  pressure  on  a spherical  array  arising  from  the  excitation 
of  a single  circular  element  has  been  obtained.  The  cases  of  a Heaviside 
step  and  a continuous  wave  (CW)  pulse  were  considered. 

The  short  time  transient  effects  were  isolated  into  efficient  creeping- 
wave  series  (Equations  28  and  1*5)  while  the  long  time  results  were  obtained 
via  classical  wave  harmonic  analyses  (Equations  9 and  10).  In  the  case  of 
a CW  pulse , the  short  time  transient  effects  were  first  separated  from  the 
steady-state  response.  The  total  solution  thus  is  represented  by  the  sum 
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of  Equations  37  and  45  or  38  and  46 


In  the  immediate  vicinity  of  the  active  element,  the  Kirchhoff  approxi 


mation  has  "been  shown  to  he  applicable.  Thus  in  this  region  the  analyses 


the  evaluation  of  the  pressure  in  the  relevant  space  and  time  domains 
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